Is there a reentrant glass in binary mixtures? 
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By employing computer simulations for a model binary mixture, we show that a reentrant glass 
transition upon adding a second component only occurs if the ratio a of the short-time mobilities 
between the glass- forming component and the additive is sufficiently small. For a ~ 1, there is no 
reentrant glass, even if the size asymmetry between the two components is large, in accordance with 
two-component mode coupling theory. For a <C 1, on the other hand, the reentrant glass is observed 
and reproduced only by an effective one-component mode coupling theory. 
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The ability to manipulate the rheological behavior of 
physical systems by the addition of suitable smaller com- 
ponents is a task of high technological interest. At the 
same time, the glassification and melting of dynami- 
cally arrested states present the theorists with a chal- 
lenging problem. Recent experiments 0,0] have revealed 
that the addition of non-adsorbing polymers to sterically- 
stabilized colloidal dispersions induces a new glass transi- 
tion scenario. Upon adding polymers to the pure colloidal 
glass, the dynamics of the colloids is first speeding up, un- 
til the glass is molten and a fluid state is reached. By 
further increasing the polymer concentration, the fluid 
vitrifles again and a reentrant glass transition material- 
izes. This new kinetically arrested phase is called "at- 
tractive" glass (in contrast to the ordinary "repulsive" 
or hard-sphere like colloidal glass), since it is generated 
by effective depletion attractions between the colloids, 
mediated by the polymers. The occurrence of an "at- 
tractive" glass seems to be generic for any system with a 
short-ranged attraction and was recently also found for 
colloids in a solvent of varying quality |3] and in micelle- 
polymer mixtures [i^. The reentrant scenario of the at- 
tractive glass had been predicted by mode-coupling the- 
ory of the colloids by treating only the colloids explicitly, 
and assuming that they interact by means of their effec- 
tive, polymer-mediated depletion interactions 5, 6, 7]. It 
has also been observed in computer simulations of parti- 
cles interacting via a short-ranged attraction |^ |9| , pro- 
vided that the range of the attraction is small enough 
with respect to the core size. Nevertheless, a general 
understanding about the mechanisms and the circum- 
stances under which additives affect the glass transition 
is still lacking. A similar situation occurs for asymmetric 
molecular mixtures (such as molten salts), but just on a 
different (microscopic) length scale. 

In this Letter, we study systematically the question of 
whether and how the dynamical properties of the addi- 
tives influence the occurrence of the glass transition, in 
the framework of a simple binary model mixture. We 
show thereby that the relative short-time mobility of the 
constituent components plays a crucial role in determin- 
ing the occurrence of vitrification of a fluid or melting of 



a glass. Our central result is that a reentrant glass sce- 
nario requires not only a large size asymmetry but also 
that the short-time mobility of the added component be 
much higher than that of the glass-forming species. In 
the case of a very high additive-mobility as compared 
to the mobility of the glass-forming component, one can 
think in terms of the adiabatic approximation: then, 
for an instantaneous configuration of the glass-forming 
component, the additives are equilibrated, establishing 
thereby the static effective interaction potential between 
the former [T3|. It is precisely in this limit that the con- 
cept of a static, effective, one-component description with 
a depletion attraction is expected to be applicable also in 
the dynamical sense. On the other hand, if the short-time 
mobilities of the two species are comparable, the deple- 
tion potential still exactly determines the partial static 
structure of the larger component j^, but not its dy- 
namical behavior. It is therefore interesting to examine 
whether the lack of validity of the depletion picture in the 
dynamical sense has an influence on the glass transition 
scenario, which is a genuinely kinetic arrest. 

From its basic formulation, the two-component mode 
coupHng theory (MOT) for the ideal glass transition 
asserts that the latter only depends on the static par- 
tial structure factors of the mixture and hence it is in- 
dependent on the individual short-time mobilities. This 
assertion holds both for Brownian dynamics (relevant for 
colloid/polymer mixtures) and for Newtonian short-time 
dynamics (relevant for molecular glass formers). Our 
computer simulation studies reveal, however, that the 
scenario and the location of the glass transition in the 
mixture depends crucially on the ratio a between the 
short-time mobilities of the glass forming component and 
of the additive. We also show that MCT correctly pre- 
dicts the glassification scenarios for the two extreme lim- 
its: if a ~ 1, the two-component MCT yields no reen- 
trance, while for a ^ 1, the effective one-component 
MCT predicts the occurrence of a reentrant glass. In 
what follows, we first describe the model and the simu- 
lation techniques. Our simulation results arc then com- 
pared to two-component and one-component MCT. Fi- 
nally, estimates for the crossover towards adiabaticity 
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and experimental consequences are presented. 

Our results are obtained for the Asakura-Oosawa (AO) 
binary model [TsL IT^ , which is a simple prototype for a 
colloid-polymer mixture but can also be used to model 
asymmetric non-additive hard-sphere mixture for molec- 
ular systems. In this model, the glass-forming particles 
("colloids") are hard spheres of diameter (Jc and the ad- 
ditives ( "polymers" ) do not interact with each other but 
only feel the colloids as hard spheres with an interaction 
radius (ctc -I- Cp)/2, where Up denotes the diameter of gy- 
ration of the polymer coils. The motivation to study such 
an asymmetric non-additive mixture is twofold: first, the 
static structure and the equilibrium phase behavior are 
known, facilitating to a great extent the understanding of 
the dynamical behavior. Second, the Asakura-Oosawa is 
a minimal model for depletion and fluid-fluid phase sep- 
aration in mixtures, thus it contains the basic physical 
mechanisms driving the occurrence of both repulsive and 
attractive glasses. Since there is no energy scale in the 
interactions, the thermal energy ksT scales out and is 
irrelevant. The remaining parameters for the AO-model 
are the size asymmetry q = Up/uc, which corresponds to 
the range of the attractive potential in the effective one- 
component picture ^3 1 ^'^d the partial number densities 
Pc Pp of the colloids and polymers, which are most conve- 
niently expressed in terms of the partial colloid and poly- 
mer packing fractions 77c = irpcCr^/Q and rjp = TTppa^/G, 
respectively. 

We carried out extensive molecular dynamics (MD) 
simulations of the AO-mixture, introducing the two 
masses nic and rrip for the colloids and the polymers. 
The latter determine the ratio of short-time mobilities 
(or thermal velocities) a between polymers and colloids 
as a = y^mp/rric- MD simulations were performed us- 
ing an event-driven algorithm [Tsj l . The relative numbers 
of particles and the volume of the simulation box were 
chosen in order to achieve equilibrium conditions for all 
studied values of q, a and r]p, within reasonable com- 
puter time jl^ . Thus, we fixed the number of colloidal 
particles to Nc = 500 for q = 0.15 and q — 0.50, and 
to Nc = 40 for q = 0.04 Partial mean squared 

displacements for colloids and polymers were calculated 
in order to estimate, via the Einstein relation, the self- 
diffusion coefficients Dc and Dp of the colloids and poly- 
mers, respectively. To reduce statistical errors, averages 
were performed over 50 independent configurations for 
q = 0.50 and q = 0.15, and over 10 for q = 0.04. The 
units of mass, length, and energy are chosen as TOc, <7c 
and fcsT. 

Simulation data for Dc as a function of the mass ratio 
are presented in Fig. ^a) for different q at a fixed state 
point. It can be seen that the diffusive motion of the 
colloids depends significantly on the mass ratio and on 
q. For sufficiently small mass ratios, the added polymer 
accelerates the dynamics if q is small and slows it down 
if q is large. This behavior is in agreement with previ- 
ous studies of the dependence of the glass transition of 
the colloids on the range of the effective attractive poten- 



9 -2 
10 



10 



10 

o 

9 
^a. 10 

9 

10' 



ri^ = 0.49 
ri' = 0.0024 



(a) 




_l I 



J I 



10' 





•-• q = 0.04 




■-■^ = 0.15 




A-A q = 0.50 






; (b) 




1 





10" 



10"' 



10 



10 



a 



FIG. 1: The dependence of the diffusion coefficients on the 
ratio of the short-time mobilities for (a) the colloidal particles 
and (b) the polymers, at the state-point rjc = 0.49, rjp = 
0.0024 and for three different g-values, normalized with Do = 
y^kBTai/rric. The horizontal line in (a) indicates the value 
of Dc when rjp — 0. Note that the legends on (a) and (b) 
apply to both panels. 



tial However, by increasing the mass ratio, the 

picture completely changes. The added polymer always 
slows down the dynamics of the colloids, irrespective of 
the range of the attractive potential. This is clear evi- 
dence that, despite of the fact that the statics is identical 
along each curve at fixed q, the dynamics is significantly 
affected by the relative mobilities of the two components. 

Data for Dp, at the same state point and varying q 
and a, are shown in Fig. ^b). For small q, Dp oc a~^, 
which pertains to the motion of an ideal gas inside the 
disordered medium formed by the colloids. For large q, 
the polymers are caged by the colloids. This causes a 
significant reduction of Dp for fixed a and a much weaker 
dependence of this quantity on a. 

In order to obtain a broader picture of the dynamical 
behavior of the system on the entire 77c — 77p-plane, we es- 
timate the isodiffusivity lines from the MD simulations. 
These can be considered as precursors of the incipient 
glass transition line and provide an estimate of the shape 
of the latter Results for the isodiffusivity lines for 

q — 0.15 are shown in Fig. [21 and for q ~ 0.50 in Fig. 
|31 In each Figure, results for two very different mass ra- 
tios, corresponding to a = 1 and a = 0.01, are shown. 
For q = 0.15, two different kinds of behavior result: for 
a = 0.01, a clear reentrance is observed, while for a = 1, 
the addition of polymers always slows down the dynam- 
ics. For a smaller size asymmetry, q = 0.50, there is 
always polymer-induced vitrification and no reentrance, 
but the impact of added polymer on the glass formation, 
embodied in the slope of the isodiffusivity curves on the 
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FIG. 2: Isodiffusivity lines on the 77c — ?jp-plane for a = 0.01 
and a = 1, accompanied by ideal glass transition lines from 
the effective one-component and the two-component versions 
of MCT. The results refer to g = 0.15. 

Vc — f^p-plane, depends strongly on a. For small mobil- 
ity ratios, the isodiffusivity curves are almost vertical, 
pointing to an insensitivity of the diffusion coefhcient on 
the amount of added polymer. For larger a, addition of 
polymer considerably slows down the dynamics. 

In order to provide a criterion for the applicability of 
MCT to different regimes of mobility ratios, we have cal- 
culated the location of the ideal glass transition within 
the framework of this theory. We have solved the im- 
plicit equations for the matrix of partial non-ergodicity 
parameters F{q), 

F{q)^S{q)~{S-\q)+nFm-' (1) 

where S{q) is the matrix of partial static structure fac- 
tors, and ^[F((7)] is the long-time limit memory ker- 
nel within MCT approximation More specifically, 
it is a quadratic functional of F(g), coupling all differ- 
ent wave-vectors, by means of coupling coefficients that 
only depend on partial static structure factors and partial 
concentrations of the system. A glass transition occurs 
when, upon changing the control parameters, the solu- 
tion of Eq. jumps from zero (ergodic state) to a finite 
value (glassy state). Note that, within MCT, inertia pa- 
rameters drop out, hence the location of the ideal glass 
line is independent of a. Eq. holds also for the one- 
component case, where S{q) and F{q) are 1x1 matrices. 

The partial structure factors needed as input for the 
two-component and one-component mode coupling the- 
ory were obtained from the recently developed fundamen- 
tal measure density functional theory [2y|- This theory 
yields analytical expressions for the partial structure fac- 
tors Sij{k) with i,j = c,p which compare well with sim- 
ulation data. In the case of the effective one-component 
MCT, S{q) = Scciq). 

Results from both effective one- and two-component 
MCT are shown in Figs. H and 01 for q = 0.15 and q = 
0.50, respectively. In the first case, the two-component 
MCT predicts a vitrification of the fiuid upon addition of 
polymers, while the effective one-component MCT pre- 
dicts a reentrant scenario. In the second case, both the 
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FIG. 3: Same as Fig.Hbut for q = 0.50. 

two component and the effective one-component MCT 
predict a vitrification of the fluid upon addition of poly- 
mers. In this last case, the slope of the ideal glass line is 
almost vertical for the one-component calculation. 

Since the isodiffusivity lines are precursors of the shape 
in the r]c — rjp plane of the glass transition line, the com- 
parison between these lines and the MCT calculations 
allows us to conclude that the MD results for small a are 
properly captured by the effective one-component MCT, 
while the MD results for large a are captured by the full 
two-component MCT calculation. 

The ratio a needed to reach the adiabatic limit-regime, 
in which an effective one-component description applies 
dynamically can be estimated as follows. During the time 
necessary for a colloidal particle to move along a distance 
CTp, there has to be a sufficiently large number of col- 
lisions, exceeding at least some typical threshold value 
A^o ~ 100. This is necessary in order to provide en ough 
statistics to feel the effective interaction dynamically [2ll . 
The resulting adiabaticity condition reads as 

a < l2Vp/{q^No). (2) 

For the two chosen two cases, a = 0.01 and a = 1 in Figs. 
121 and 13 we realize that, for a typical rjp, the abiabatic- 
ity criterion Q is fulfilled in the former but not in the 
latter case. This in turn explains why the effective de- 
pletion potential is dynamically meaningful for the case 
a — 0.01, so that the effective one-component MCT is 
applicable. On the other hand, for a = 1 adiabaticity 
is not achieved and both the colloids and polymers are 
cage formers, so that the two-component MCT has to be 
employed in the study of the vitrification transition. 

The same scaling arguments employed for the deriva- 
tion of criterion jSJ can be carried out for short-time 
Brownian dynamics, for which a = / with ui"'' 
and Dp"'' denoting the short-time diffusion coefficients 
of the colloids and the polymers, respectively, in their 
common solvent. In a mixture of uncharged colloids and 
polymers, these coefficients scale with the inverse radii 
of the particles, according to Stokes' expression, hence 
they are coupled to q via d'^'^ /D^^ — q. Then the 
critical asymmetry q below which adiabaticity holds is 
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given by Qc = (12r7p/7Vo)^/'^, which is close to 0.1 for 
typical values of rjp and Nq, On the other hand, in mix- 
tures of charged suspensions, the physical hard core di- 
ameters ai^^ and dp"' are different and independent from 
the effective interaction diameters that enter in the q ra- 
tio. Indeed, effective interaction diameters (Tc and ap are 
dictated by the long-range Coulomb interactions, thus 
they can greatly exceed the physical ones. The charge 
on the particles provides in this case the physical pa- 
rameter that allows tuning of the physically relevant size 
ratio q. In this case, the adiabaticity criterion Q reads 

as 4°V'^c"' < l27jp/iq^No). 

In concluding, we have shown that a sufficient asym- 
metry in the short-time mobilities is necessary for the 
reentrant glass-scenario to materialize in highly asym- 
metric binary mixtures. For glass- forming mixtures gov- 
erned by molecular dynamics, the glass transition de- 
pends on the mass ratio between the two components 
[23, but huge mass asymmetries are needed to achieve 
the adiabatic limit. Vastly different mass ratios can be 
realized in binary mixtures of dusty plasmas whose 
dynamics is almost Newtonian. On the basis of the con- 
siderations put forward in this Letter, the size asymme- 
try q for colloid-polymer mixtures has to be smaller than 
0.1 in order to get a reentrant glass. This condition is 
stronger than the one obtained by using a priori an effec- 
tive one-component picture. Indeed, in the present cal- 
culation one-component MCT would predict a reentrant 
glass for q < 0.3, while for previous studies of effective 
one-component pictures 0, [3 q up to at least 0.2 was 



found to have an attractive glass |2J|. This finding helps 
to explain why in experimental studies a very small 

size asymmetry was indeed necessary in order to obtain 
an attractive glass. The most striking consequences are 
obtained for a mixture of charged colloids: a reentrant 
glass can be lost if the physical core of the high-charge 
particles is decreased, even if the charges are kept con- 
stant. By tuning the short-time mobility, the glass for- 
mation upon addition of a second repulsive component 
can be tailored, opening a way for external manipulation 
of the rheology of the mixture. 

If the kinetic glass transition in a mixture is calcu- 
lated within mode coupling theory, it is not known a 
priori whether the two-component or the effective one- 
component versions of the theory has to be employed. 
We have shown that it is the adiabaticity criterion (0) 
that determines which theory must be applied. The im- 
pact of the a-dependence of the glass transition should 
also have unexplored important consequences in polymer 
mixtures, such as, e.g., mixtures of star polymers of dif- 
ferent arm numbers and radii, in which the lighter poly- 
meric component has been used as a mo dify ing agent of 
the flow-behavior of the heavier one [23, f23| . 
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